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We apply high-order many -body perturbation theory for the calculation of ground-state energies of closed-shell nuclei using realistic 
nuclear interactions. Using a simple recursive formulation, we compute the perturbative energy contributions up to 30th order and 
compare to exact no-core shell model calculations for the same model space and Hamiltonian. Generally, finite partial sums of 
this perturbation series do not show convergence with increasing order, but tend to diverge exponentially. Nevertheless, through a 
simple resummation via Pade approximants it is possible to extract rapidly converging and highly accurate results for the ground 
state energy. 
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1. Introduction 

\ The treatment of the nuclear many-body problem is a cen- 
tral and long-standing issue in nuclear structure theory. Ideally, 
we would like to solve the many-body problem ab initio, i.e., 
starting from a given nuclear Hamiltonian without any concep- 
tual approximations. With the advent of high-precision nuclear 
J potentials that are based systematically on Quantum Chromo- 
dynamics (QCD) through chiral effective field theory Hid me 
demand for exact ab initio solutions of the nuclear many-body 
problem has grown. Only these schemes establish a rigorous 
and quantitative connection between nuclear structure observ- 
ables and the underlying QCD input. 

The no-core shell model (NCSM) is one of the most univer- 
sal exact ab initio methods, which gives access to all aspects of 
nuclear structure OHHQ]. Other methods, are either restricted 
to certain classes of Hamiltonians, like the Green's Function 
Monte Carlo approach @|, or they are limited to certain nuclei 
and observables, like the coupled-cluster approach (Q]. All of 
them are computationally demanding, which leads to a severe 
limitation regarding the number of nucleons that can be han- 
dled. 

Therefore, approximate many-body schemes using the same 
Hamiltonians, i.e. approximate ab initio methods, also pro- 
vide indispensable information. In particular approaches that 
use controlled and systematically improvable approximations 
are of great practical importance. In this category, many-body 
perturbation theory (MBPT) is one of the most powerful and 
widely used methods. On the one hand, the evaluation of low- 
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orders of perturbation theory is computationally simple and can 
be done for the whole nuclear mass range fl H TToL U ] as 
well as for infinite nuclear matter ILL 



is deemed systematically improvable, either by extending the 
MBPT calculations order-by-order or by using infinite partial 
summations, like ladder- or ring-type summations lfl3i Ua. [l5ll . 
However, the accuracy of low-order perturbative estimates, e.g. 
for ground-state energies, or possible extensions of the MBPT 
series to higher orders and the resulting convergence pattern are 
rarely, if ever, addressed in the nuclear structure context. 

In this paper, we apply MBPT for the calculation of the 
ground state energy of several closed-shell nuclei. We extend 
the order-by-order calculation of the perturbative energy contri- 
butions up to 30th order, study the convergence behavior, and 
compare to exact NCSM calculations for the same Hamiltonian 
and model space. We introduce Pade approximants as a highly 
efficient tool for the resummation of the divergent power-series 
of MBPT into a rapidly converging series and demonstrate their 
accuracy for the description of ground-state energies. 

2. Many-body perturbation theory 

2.1. Formalism 

We aim at a perturbative expansion of the many-nucleon 
Schrodinger equation 



HJ¥n)=E n PF„> 



(1) 
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for the translational invariant nuclear Hamiltonian H — T — 
Tcm + V, where we assume V to be a two-body interaction for 
simplicity. In a first step we have to chose the unperturbed ba- 
sis, which in turn defines the unperturbed Hamiltonian. From 
the practical point of view, a basis of Slater-determinants con- 
structed from a set of single-particle states is most convenient. 
The underlying single-particle basis will typically be a Hartree- 
Fock or a harmonic oscillator basis — for simplicity we assume 
the latter. The unperturbed Hamiltonian Hq is a one-body oper- 
ator containing the kinetic energy T and a harmonic oscillator 
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potential. The unperturbed Slater determinants |0„) fulfill the 
eigenvalue relation 



which characterize the perturbative corrections to the eigen- 
states l^f') expanded in the unperturbed basis 



H \®„) = e n |O n ) 



(2) 



with eigenvalues e„ being the sum of the single-particle ener- 
gies of the occupied states. After the unperturbed Hamiltonian 
is fixed, the perturbation is defined through W — H — //(>. This 
partitioning leads to the M0ller-Plesset formulation of MBPT 
and obviously other partitionings of the Hamiltonian are pos- 



sible 11161 11711 . For ease of presentation, we assume that the 



unperturbed state corresponding to the eigenstate we are inter- 
ested in is non-degenerate, as it is the case for the ground state 
of closed shell nuclei. In the case of degeneracy, as e.g. for 
the excited states of closed shell nuclei, one would have to di- 
agonalize the full Hamiltonian in the degenerate subspace and 
pick the eigenstates with the desired quantum numbers as un- 
perturbed states. 

The standard Rayleigh-Schrodinger perturbation series can 
now be constructed based on a Hamiltonian (using the notation 
from Ref. 0) 

H(A) =H + XW (3) 

containing an auxiliary expansion parameter A that continu- 
ously connects the unperturbed Hamiltonian Hq = H(A = 0) 
with the full Hamiltonian H = H(A =1). The energy eigenval- 
ues E„(A) and the corresponding eigenvectors |*P„(/l)) of H{A) 
are formulated as a power series in A 



E n (A)=E^+AE^ + A^> + 



\¥ n (A))= \yf>) + A\V n l >) + A'\V%>) + 



(4) 



In the absence of degeneracy the lowest-order contributions are 
simply given by the unperturbed quantities, i.e., 

E ( °> = e„ , yVf>) = |0„) . (5) 

Inserting the Hamiltonian (0) and the power series (0| into 
the Schrodinger equation (HJ leads to the fundamental equation 



#0 |¥f > + J] A "( W ^n'^) + Ho 

oo p 



(6) 
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Assuming that the unperturbed states form an orthonormal basis 

■AO) 



1^) 



(9) 



with C%1 = for p > and = 5 njn . 

We can cast Eqs. (|7]i and ([8) into a more transparent form by 
systematically introducing the amplitudes C { „ p ^ n and formulating 
all matrix elements in terms of the unperturbed states. For the 
pth-order energy contribution we obtain 

=Yj&n\W\® m )C%- l) . (10) 

m 

Similarly we obtain for the pth-order amplitudes 



C 



(p) 



-(X<<l> m |W|cD ^l/ )C, ( ^ 1, 
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Together with C, ( , ( j„ = and E ( ^ = e„ these relations form a 
recursive set of equations which uniquely determines the per- 
turbative corrections for all energies and states to all orders. 

Usually one would use these general expressions to derive 
explicit formulae for the lowest-order corrections. The ma- 
trix elements of the perturbation in the unperturbed Slater- 
determinant states can be evaluated explicitly and the summa- 
tions over the many-body basis set can be replaced by summa- 
tions over single-particle states. In this way we would recover 
the standard expressions for, e.g., the second- and third-order 
energy corrections liifTllflal. 

2.2. Evaluation to high orders 

When attempting to evaluate the perturbative corrections be- 
yond third- or forth-order the explicit formulae for the energy 
corrections become impractical because of the large number of 
nested summations. A much more elegant way to evaluate high- 
order contributions makes use of the recursive structure of Eqs. 
( TTOl ) and (fTTT) . The only ingredients needed are the many-body 
matrix elements of the full Hamiltonian H with respect to the 
unperturbed basis |<t>„). Starting from the zeroth-order coeffi- 



cients C 



((» 



Snjn we can readily evaluate the first-order en- 



ergy contribution E^ from ( [Tol l. This in turn allows us to com- 
pute the first-order coefficients C„^ via (fTTT i. Generally, for the 



and using the intermediate normalization (% '\*¥ n {A)) = 1 we evaluation of the energy contribution E„ only the coefficients 



obtain QVT^f) = for p > 0, which allows us to project-out c% n l) of the previous order are required. For the evaluation of 



all required information on the individual contributions in the 
power series. By multiplying Eq. (O with (4^ '| and matching 
same orders of A on both sides, we immediately obtain a simple 
expression for the pth-order energy contribution 



E^ = (XI . 

By multiplying Eq. © with (\Pj£ 



(7) 



with m + n and matching 
/l-orders, we obtain an expression for the amplitudes 



C 



7(0) 
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7=1 



(0)|xp(p-7) 



>) 



(8) 



the coefficients C { ,j'm all energy contributions up to order p and 
all coefficients up to order (p - 1) need to be known. 

Technically, the recursive evaluation of the perturbation se- 
ries bears some resemblance to the Lanczos algorithm for the 
iterative solution of the eigenvalue problem for a few extremal 
eigenvalues as it is used in the NCSM. The most significant op- 
eration is a matrix- vector multiplication of the Hamiltonian ma- 
trix with the coefficient vector from the previous order, which 
constitutes the first term in the evaluation of the coefficients 
(fTTT) . Because the second term in (fTTI) involves the coefficient 
vectors from all previous orders, we store them for simplic- 
ity. These computational elements are the same as for a simple 



Lanczos algorithm in the NCSM or in corresponding configura- 
tion interaction (CI) approaches, therefore, an implementation 
of high-order MBPT using NCSM or CI technologies is straight 
forward. However, since the computational elements are the 
same, so are the computational limitations: This direct imple- 
mentation of high-order MBPT is limited to the same model 
spaces and particle numbers as the full NCSM. This is not a 
concern for the present study, but for an application of MBPT 
beyond the domain of the NCSM one has to resort to other eval- 
uation schemes. 

2.3. Applications: A He, 16 <9 and 40 Ca 

As examples for a direct application of high-order MBPT 
we consider the ground-state energies of 4 He, 16 0, and 40 Ca. 
Throughout this work we use an intrinsic Hamiltonian with 
a soft two-nucleon interaction that is derived from the chiral 
N3LO potential (H via a Similarity Renormalization Group 
(SRG) transformation OH EH. The final flow parame- 
ter for the SRG evolution of the interaction is a — 0.02 fm 4 
which corresponds to a momentum scale of A = 2.66 fm . 
This choice for the flow parameters leads to a unitarily trans- 
formed interaction which is sufficiently soft to warrant excel- 
lent convergence properties with respect to model space size in 
the NCSM but at the same time yields ground-state energies 
which are in reasonable correspondence with experiment in the 
mass-range considered here. To allow for a direct comparison 
with exact NCSM calculations for the same Hamiltonian and 
the same model space, we use an A^ max ftf2 model space also for 
the MBPT calculations. We have confirmed, however, that all 
conclusions regarding the performance and limitations of the 
MBPT approach do not depend on this particular choice. 

In Fig. Q] we summarize the results of an order-by-order 
MBPT calculation up to 30th order for the ground state energy 
of the three nuclei. For 4 He the calculations were performed in 
a 12 fill model space, for lfi O in 6HQ., and for 40 Ca in AMI, each 
with two different oscillator frequencies Ml. The partial sum of 
the perturbative energy contributions up to order p, 



p 

p'=0 



(/>'> 



(12) 



is depicted in upper row and the modulus of the individual pth- 
order contributions, |£^|, on a logarithmic scale in the lower 
row. Here and in the following we omit the index n = for 
convenience. 

Already the first glance at Fig.[T]reveals a fundamental prob- 
lem with the convergence behavior of the perturbation series. 
For 4 He, as depicted in Fig. [Ha), we observe two different pat- 
terns depending on the oscillator frequency. For Ml = 20 MeV 
the partial sum E sma (p) shows an alternating behavior with a 
systematically decreasing amplitude. Beyond 10th order one 
might consider the perturbation series converged and the result- 
ing energy is in excellent agreement with the result of an NCSM 
calculation with the same Hamiltonian in the same model space. 
However, a change of the oscillator frequency destroys this pic- 
ture. For Ml = 32 MeV, where the NCSM provides a lower 
ground-state energy, we again observe an alternating sequence 



of energy contributions E (p \ but this time without any sign of 
convergence. The absolute value of the individual energy cor- 
rections does not decrease with increasing order, it even shows 
a slightly increasing trend. Hence even in the simplest case, the 
4 He ground state, the convergence of the perturbation series is 
not guaranteed. 

The situation is even more dramatic for lfi O or 40 Ca as de- 
picted in Fig.QIb) and (c). In all cases the size of the perturba- 
tive energy contributions grows exponentially with p. The 
partial sum E sam (p) exhibits a strong oscillatory behavior with 
increasing amplitude. Only the lower orders, typically up to 
10th order for 16 and up to 5th order for 40 Ca, lead to binding 
energies in a physically meaningful energy range. At the 30th 
order the perturbative contributions are in the order of 10 4 MeV 
for lfi O and 10 9 MeV for 40 Ca — this is beyond any physical en- 
ergy scale present in the nuclear many-body problem. We were 
not able to find a convergent scenario by varying the oscillator 
frequency or the model space size or truncation for these nuclei. 

The explosion of the perturbative corrections beyond any 
meaningful energy scale suggests a principal mathematical 
problem in the representation of the energy eigenvalue E(A) as 
a partial sum of a simple power series (|4]i . 

3. Pade approximants 

3.1. Formalism 

Prompted by the drastic failure of a partial sum of a simple 
power series to describe the energy E(A) at the physical point 
A — 1 one might consider more general expansions of this func- 
tion. A next step would be an expansion of the energy E(A) in 
terms of a rational function composed of separate power series 
for numerator and denominator 



E(A) 



A(A) a + a\A + a 2 A + 



B(A) b + b 1 A + b 2 A 2 + 



(13) 



Obviously we will not attempt to re-derive perturbation theory 
for this type of expansion. The above is useful only, if we could 
use the information contained in the standard MBPT energy 
contributions E (p) to construct this rational expansion. 

Exactly this is achieved through the Pade approximants 121 
2211 . Given a power series (0]i of the function E(A), then the 
Pade approximant 



M 



[M/N](A) = 



«o + o\A + a 2 A + ■ • ■ + omA 
b +b l A + b 2 A 2 + --- + b N A N 



(14) 



with numerator being a polynomial of order M and the denomi- 
nator a polynomial of order N is constructed such that its Taylor 
expansion reproduces the first M + N orders of the initial power 
series, i.e. 

E{A) = [M/N](A) + 0(A M+N+1 ) . (15) 



From this definition one can immediately construct a coupled 
system of equations that determines the coefficients a„ and 
b m of the Pade approximant from a given set of E (p) with 
p — Q,...,N + M. An alternative and more elegant form 



1 2JJ, |22|] relates the Pade approximants to determinants of two 
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(a) 4 He 



(c) 40 Ca 




Figure 1: Contributions to the ground-state energies in MBPT up to 30th order for different nuclei and model spaces: (a) He in a Y2M1 model space with Ml = 20 
MeV (•) and 32 MeV (♦); (b) lfi O in 6M1 with Ml = 20 MeV (•) and 24 MeV (♦); (c) 40 Ca in 4M1 with Ml = 20 MeV (•) and 24 MeV (♦). The upper 
panels depict the partial sum E sum (p) defined in Eq. )12t as function of the largest order p, the lower panels the modulus of the individual contributions £' p ' on a 
logarithmic scale. The dashed horizontal lines in the upper panels indicate the NCSM ground-state energies for the respective nuclei and model spaces. 



(N + l)x(N + 1) matrices containing directly the power-series 
coefficients E (p) 



[M/N]W : 



£(M-N+1) 


£(M-N+2) 


£(M+1) 


£(M-N+2) 


£(M-N+3) 


£(M+2) 




£(*»+!) 


E (M+N) 


M-N 


M-N+l 


M 


2 E®A N+i 


2 E {,) A N+, - [ ■ 


■ ■ 2 E®X 


i=0 


i=0 


1=0 



£(M-N+l) £(M-N+2) 
£(M-N+2) £(M-N+3) 



E (M) 
A N 



A N-l 



£(M+l) 

£(M+2) 



£{M+N) 
1 



(16) 

where we set Zs (p) = for p < 0. We will use this form to 
evaluate various Pade approximants in the following. 

Before considering numerical results, we should like to men- 
tion a few formal properties of the Pade approximants that 
are of importance for the present application. The mathemat- 
ical foundation for using Pade approximants for our purpose 
in the first place is provided by the Pade conjecture (simpli- 
fied) lull 12211 : Let E(A) be a continuous function for \A\ < 1, 
then there is an infinite subsequence of diagonal Pade approx- 
imants [N/N](A) that for N — > 00 converges locally uniformly 
against E(A) for \A\ < I. For our application the continuity re- 
quirements for the function E(A) are always fulfilled, thus we 
expect the diagonal Pade approximants to show a convergence 
behavior — unlike the simple power series. 

Additionally the Pade approximants have a number of spe- 
cific properties that would be extremely valuable in the present 
context. If the power series expansion of E(A) is a Stieltjes se- 
ries, then the Pade approximants fulfill the condition 



for A > as well as a whole set of related inequalities 12 ll 12211 . 
Thus the diagonal and the super-diagonal Pade approximants 
provide upper and lower bounds for the full energy E{A), re- 
spectively. Furthermore, these bounds improve monotonically 
with increasing order M of the Pade approximant. Unfortu- 
nately, it turns out that the power series we start from is not a 
Stieltjes series in general. 

3.2. Applications: A He, l6 0and 4a Ca 

Using the results of the order-by-order calculation of the en- 
ergy corrections E (p ^ up to 30th order of MBPT we can con- 
struct all Pade approximants with N + M < 30 from Eq. (1161 l. 
Evaluating the approximant at A = 1 yields an estimate for the 
ground-state energy of the perturbed system 



£pade(M/A0 = [M/N](A = 1) 



(18) 



[M/M](A) > E(A) >[M- l/M](A) 



(17) 



We will focus on the diagonal Pade approximant, Ep. d fe(M/M), 
and the super- and sub-diagonal approximants, Ep^eiM - 1 /M) 
and £'p a de(M/M - 1), respectively, because of the convergence 
and boundary theorems available for those. 

A collection of all diagonal as well as sub- and super- 
diagonal approximants with N + M < 30 for 4 He, le O, and 40 Ca 
using the oscillator frequencies that yield the lowest ground- 
state energy is provided in Fig. [2] The first remarkable ob- 
servation is that the Pade approximants converge very quickly 
for sufficiently large order — we have observed this behavior in 
all cases we considered. For M + N > 10 essentially all Pade 
approximants provide the same ground-state energy. We em- 
phasize that the input for the construction of those Pade ap- 
proximants are the exponentially diverging coefficients from the 
power-series formulation of MBPT discussed in Fig. Q] The 
Pade resummation of these coefficients efficiently regularizes 
these divergencies and leads to exceptionally stable results for 
all orders M + N > 10. 
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(a) 4 He 



(b) 16 



(c) 40 Ca 
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Figure 2: Pade approximants for the ground-state energies of (a) He, (b) O, and (c) Ca as function of the summed order M + N. The different symbols represent 
the diagonal approximants Ef^(MIM) (•), the super-diagonal approximants £p a de(^ — (♦). and the sub-diagonal approximants Ep^i(M/M - 1) (A). The 
model space size N lTmx and the oscillator frequency Ml is quoted in the individual panels. The dashed horizontal lines indicate the NCSM ground-state energies for 
the respective nuclei and model spaces. 



The second important observation results from the compari- 
son of the converged Pade approximant^j] with the exact energy 
eigenvalue obtained from a solution of the matrix eigenvalue 
problem for the Hamiltonian in the same model space — i.e., 
from the corresponding NCSM calculation — as indicated by the 
dashed horizontal line in Fig. [2] The converged Pade approxi- 
mants exactly reproduce the corresponding energy eigenvalues, 
i.e., Pade resummed perturbation theory and the exact solution 
of the eigenvalue problem become equivalent. 

A quantitative comparison is presented in Tab. [1] where the 
difference of the partial sums E sum (p) and the Pade approxi- 
mants E-pzfeiM IN) to the exact NCSM eigenvalues £ncsm are 
shown. The latter were obtained using the Antoine shell-model 
code [23]. Starting from M + N » 10 the deviations of the Pade 
approximants from the exact result are getting very small and 
starting from M + N ~ 20 the Pade approximants are numeri- 
cally identical to the exact result for all nuclei. In this regime 
the individual MBPT contributions E (p) are already increasing 
exponentially for 16 and 40 Ca (cf. Fig. [TJ and the partial sum 
Esum(p) does not provide any sensible estimate of the ground- 
state energy. The Pade approximants prove to be a highly ef- 
ficient tool to extract a virtually exact and stable result for the 
energy from the first 10 or more coefficients E (p) of the strongly 
fluctuating and non-converging power series. Considering the 
scale of the order-to-order fluctuations and the absolute size of 
the perturbative contributions E (p) the stability and the precision 
of the converged Pade approximants is truly remarkable. 

For application purposes, the behavior at low orders is also 
of interest. As shown in Tab. [T| the deviations AE sum (p) and 
&Ep ai ie(M/N) are of comparable magnitude up to about fifth or- 
der, both showing sizable fluctuations. Hence, in this low-order 
domain the Pade approximants do not improve on the results 
obtained from a simple partial sum. Only beyond M + N ~ 5 
do the Pade approximants start to converge, i.e., the variations 



Here, the term convergence refers solely to the convergence with respect 
to the order M + N and not to convergence with respect to the model-space size 
Mnaxi which is a separate issue. 



within a set of approximants of neighboring order reduce sys- 
tematically. At the same time the deviations of the partial sums, 
AE sam (p), start to increase exponentially. 

The stability of the Pade approximants Ep. d fe(M/N) across 
various neighboring orders M and N is an important intrin- 
sic criterion for convergence and for the accuracy of the Pade 
approximants as compared to the exact result. Therefore, it 
seems advisable to always consider sets of several approxi- 
mants. Moreover, there is always the possibility that individual 
approximants completely escape the general trend, such as the 
£pade(2/2) approximant for 40 Ca at Ml = 20 MeV that has a 
large positive and thus unphysical value. These cases are a re- 
minder that the convergence theorems for Pade approximants, 
e.g., the Pade conjecture, only cover subsequences of approx- 
imants. Finally we note that the MBPT power series in the 
present examples turns out not to be a Stieltjes series. As the 
Pade approximants of Tab.Q]show, the inequality ( fTTl i as well as 
related inequalities are not fulfilled. The Pade approximants for 
the ground state energy on nuclei in the present MBPT frame- 
work do not provide rigorous bounds for the exact eigenvalues. 

4. Conclusions & Outlook 

We have formulated and applied many-body perturbation 
theory up to high orders for the description of ground-state en- 
ergies of closed-shell nuclei using realistic Hamiltonians. In 
contrast to typical applications of MBPT in nuclear physics that 
are limited to second or third order, we extend the order-by- 
order evaluation of the perturbation series up to 30th order using 
a simple recursive scheme. In order to facilitate the comparison 
with exact eigenvalues obtained in NCSM calculations, we have 
limited ourselves to an harmonic-oscillator single-particle basis 
and a N mdx Ml space. However, results for other single-particle 
bases and model-space truncations are qualitatively similar. 

Our major results and conclusions are: First, a simple par- 
tial sum of the perturbative series in general does not converge. 
On the contrary, we typically observe an exponential increase 
of the individual perturbative contributions \E™\ and an oscilla- 
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Table 1 : Large-scale MBPT results for the ground-state energies of 4 He, lb O, 
and 40 Ca. Shown are the exact NCSM energies for the respective model space 
and the deviations of the partial sums, A£ sum (p) = E sum (p) - £ncsm, as well 
as the deviations of various Pade approximants, AEp^iM/N) = Ep^iM/N) - 
£ncsm ■ All energies are given in units of MeV. 
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+ 14.93 


+ 13.06 


+3115.69 


-193.45 


AFd i-(2/3) 


+0.108 


+0.892 


+8.19 


+9.19 


+29.10 


+61.98 


AFd W3/2^ 


+0.066 


+0.865 


+ 1.01 


+ 1.32 


+3.61 


+9.34 


AFd 


+0.047 


+0.761 


-1.32 


-2.34 


-6.00 


-16.04 


AFd i-n/4) 


+0.008 


+0.066 


+0.66 


+ 1.48 


+2.34 


+4.01 


A£ Pad e(4/3) 


+0.132 


+2.666 


-0.11 


-0.11 


-0.88 


-2.96 


A£ Pad i(4/4) 


+0.019 


+0.492 


+0.18 


+0.05 


+3.45 


+26.11 


A£ Pad e(4/5) 


+0.015 


+0.314 


-0.80 


+0.03 


+2.43 


+6.46 


A£pade(5/4) 


+0.012 


+0.135 


+0.03 


+0.03 


+0.22 


+ 1.14 


A£ Pad i(5/5) 


-0.029 


+0.136 


-0.04 


+0.05 


-0.50 


-2.45 


A£ Pad i(6/6) 


-0.037 


+0.040 


+0.01 


+0.00 


+0.16 


+ 1.37 


A£pade(8/8) 


+0.001 


+0.001 


+0.01 


+0.01 


+0.01 


+0.20 


AEp^^lO/lO) 


+0.000 


+0.001 


+0.01 


+0.01 


+0.00 


+0.01 


A£pad<(12/12) 


-0.000 


-0.002 


+0.01 


+0.01 


+0.00 


+0.01 


A£pade(15/15) 


-0.000 


-0.000 


+0.01 


+0.01 


+0.01 


-0.04 



tory behavior of the partial sum E sum (p) as function of p. Thus, 
finite partial summations, even if they are extended to high or- 
ders, do not provide a stable and systematically improvable ap- 
proximation for the exact energy eigenvalue. 

Second, Pade approximants offer a computationally simple 
yet powerful tool to extract a convergent series from a finite 
set of perturbative energy contributions Solely through a 
resummation of the finite /9th-order power series to a rational 
function, whose Taylor expansion up to order p is identical to 
the initial power series, we are able to extract a highly stable and 
convergent approximation for the energy. The information en- 
tering these Pade approximants of order p — M + N is identical 
to the power series of order p and so is the computational effort. 
However, whereas a finite partial sum E smn (p) explodes with in- 
creasing order p, the Pade approximants Ep a fe(M/N) converge. 

Third, beyond a sufficiently large order, typically M + N > 
10, the different Pade approximants become very stable and 
converge rapidly to a unique value for the energy. This energy 
is identical to the exact eigenvalue obtained in the correspond- 
ing NCSM calculation. In this sense, Pade resummed MBPT 
and the Lanczos diagonalization become numerically equiva- 
lent. Unfortunately, the computational effort is also comparable 
at least for the implementation we adopted for the MBPT here. 

Fourth, at low orders, i.e. M + N < 4, the Pade approx- 



imants generally do not yield an improved approximation for 
the energy. Their deviations from the exact eigenvalue fluctu- 
ate and are of the same order of magnitude as the errors of the 
corresponding partial sums. If one is limited, for computational 
reasons, to very low orders, then the version of MBPT used here 
can only provide a rough estimate that might differ significantly 
from the exact eigenvalue. 

These initial studies open a number of new avenues for the 
study and application of MBPT in nuclear structure. Beyond 
the MBPT calculations presented here, one can use optimized 
single-particle bases and improved partitionings of the Hamilto- 
nian to influence the convergence behavior of a finite order-by- 
order MBPT calculation. Furthermore, one can exploit infinite 
partial summations, e.g., ladder- or ring-type summations in- 
cluding MBPT contributions from all orders, and compare their 
quality to Pade resummed finite-order calculations. These im- 
provements and extensions will be the subject of future studies. 
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